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' We expand a partial diflFerence équation (PAE) on multiple lattices and obtain the PAE which 

governs its far field behaviour. The perturbative-reductive approach is here performed on well known 
nonlinear PAEs, both integrable and non integrable. We study the cases of the lattice modified 
»^ Korteweg-de Vries (mKdV) équation, the Hietarinta équation, the lattice Volterra-Kac-Van Moer- 

f| . beke (VKVM) équation and a non integrable lattice KdV équation. Such réductions allow us to 

ÇL(' obtain many new PAEs of the nonlinear Schrôdinger (NLS) type. 
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1 Introduction 



Problems involving the évolution of nonlinear phenomena, both continuons and discrète, have become of 
increasing interest in varions branches of science and engineering. Nonlinear waves, without dissipation 
and dispersion give rise in a finite time to a discontinuity. A typical example of nonlinear wave is the shock 
wave produced by a supersonic object. Dissipation and dispersion play an important rôle in balancing 
the steepening due to nonlinearity, so that when thèse effect are présent, a steep but smooth solitary 
wave may be formed and then propagates for ail times. The solitary wave phenomenon has actually 
been observed for many years in the form of a surface wave in shallow water. A model équation of 
nonlinear dispersive phenomena may, in gênerai, be very complicate. The soliton may appear only in 
the asymptotics, after a long transient period. Thus to be able to put in évidence the solitons, Taniuti 
and coUaborators [T^ introduced an asymptotic method which makes it possible to reduce gênerai 
nonlinear évolution équations to some more tractable nonlinear équations. This method go under the 
dénomination of reductive perturbation technique. Under the assumption that the amplitude of the waves 
are small, one is able to reduce the starting hyperbolic System to a few simple équations, such as the 
Burgers équation, the Korteweg-de Vries équation, the nonlinear Schrôdinger équation and few others. 

In the reductive perturbation method, the space and time coordinates are stretched in terms of a 
small expansion parameter and we introduce the concept oî far field, as the field governing the asymptotic 
behaviour of the reduced équation. To give a simple idea of the reasoning underlining this concept, let 
us consider, as an example, the familiar wave équation in two variables: 

4>,tt - (j>.xx = 0. (1) 

The gênerai solution of équation JQ) can be expressed as the superposition of waves moving to the right 
and to the left. In gênerai thèse two waves are excited simultaneously by an arbitrary initial condition. 
However, if the initial condition is localized, after a certain finite time the disturbance séparâtes in a 
progressive wave propagating to the right and one to the left, and they are solutions to a first order 
équation, i.e. an équation of one fewer degree of freedom: 

0,x ± 4>,t = 0. 

We call the solutions of the first order équation the far field solutions of the original wave équation. The 
concept of far field came from the idea of finding properties of a given évolution équation which do not 
dépend in a sensitive manner on the détails of the initial conditions, but correspond to a wide class of 
initial conditions. 

As an example of the simplification obtained by considering the reductive perturbation method, let 
us consider a Riemann wave: 

0,t+A(0)0,, =0. (2) 

When the wave function is small we may find the solution by a perturbation calculation. Let e be a 
small parameter and let us expand the solution around the constant solution 



Expanding in powers of e we get from équation (0) the foUowing results: 



e 



. 



with 



Ao = A(</.(°)), A,,,o, = (^^) 



0=0(0) 

Introducing the new variables x' — x — X^t, t' — et, we can rewrite the équation ((SJ, up to the second 
order in e, as: 

</>W+A,^(o,</.«05 =0. (3) 

If we consider a nonlinear dispersive System, like, for example, the Euler équation, instead of équation 
© we get: 

4>,t' = (ji^x'x'x' + Q(t>4>,x': (4) 
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that is the Korteweg-de Vries (KdV) équation. 

The nonlinear System at the lowest order approximation can admit a solution given by monocromatic 
wave packets, i.e. — Aexp[kx — Lij(k)t]. Than it is reasonable to consider perturbations of such 

solution and to turn the nonlinear System into a set of équations for the complex envelope of thèse 
packets. The characteristic packet size and wavelenght play the rôle of différent scales for this System. 

Let us consider, for example, the KdV équation Q for a small amplitude field (p of order e. The lincar 
équations admits a monocromatic solution with dispersion relation (jj{k) = —k^ . Then the solution of the 
KdV équation can be written as 

+00 

n— — oo 

with 

x' = e{x + 3kh), t' ^ -6e^kt, 
«0 = 1, a„ = n — 1, n > 1, 
and vi will satisfy the well known integrable nonlinear Schi^odinger (NLS) équation 

viA' + ^vi^x'x' - k'^vi\vi\'^ ^0. (5) 

It is important to notice that thèse multi-scale expansions are structurally strong and can be applied 
to both integrable and non integrable Systems. Zakharov and Kuznetsov in the introduction of their 
article jl4j say: // the initial System is not integrable, the resuit can be both integrable and nonintegrable. 
But if we treat the integrable System properly, we again must get from it an integrable System. 

Calogero and Eckhaus 3 used similar ideas starting from generic hyperbolic Systems to prove in 1987 
the necessary conditions for the intcgrability of nonlinear partial differential équations (PDEs). Later 
Degasperis and Procesi .4 introduced the notion of asymptotic integrability of order n by requiring that 
the multi-scale expansion be verified up to order n. 

Also in the case of differential équations on a lattice, we would like to have a reliable reductive 
perturbative method which would produce reduced discrète Systems. As the far field solution implies the 
introduction of a new variable which combines the continuons time with the discrète lattice, it is natural 
to get from a differential-difference équation by the reductive perturbation technique a continuons NLS 
équation © • Léon and Manna |Zj and later Levi and Heredero 9,^ proposed a set of tools which allows to 
perform multiscale analysis for a discrète évolution équation. Thèse tools rely on the définition of a large 
grid scale via the comparison of the magnitude of related différence operators and on the introduction 
of a slow varying condition for functions defined on the lattice. Their results, however, are not very 
promising as the reduced modcls arc ncither simpler nor as integrable as the original ones. Starting from 
an integrable modcl, like the Toda lattice, Léon and Manna [2| produce a non integrable differential- 
difference équation of the discrète NLS type. Levi and Heredero [H] from the integrable differential- 
difference NLS équation got a non integrable System of differential-difference équations of KdV type. 

In the présent paper we consider the case of nonliner partial différence équations (PAEs). To be able to 
carry ont the discrète reductive perturbation technique, in section|21we introduce multiple lattice variables 
and give a définition of slow varying functions on the lattice. Section |3|is devoted to the application of 
the perturbative expansions introduced to the case of a set of integrable and non integrable équations, i.e. 
the lattice modified Korteweg-de Vries (mKdV) équation, the Hietarinta équation, the lattice Volterra- 
Kac-Van Moerbeke (VKVM) équation and a non integrable lattice KdV équation. Section 0]is devoted 
to some conclusive remarks. 

2 Multiple— scales on the lattice and functional variation on them 

The aim of this section is to fix the notation and to introduce the mathematical formulae necessary to 
reduce integrable and non integrable lattice équations in the framework of the perturbative-reductive 
approach. In doing so we will partly foUow [H], trying to présent a clearer and simpler dérivation of ail 
necessary formulae. 
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2.1 Slow varying variables on the lattice 

Given a lattice defined by a constant lattice spacing h, we will dénote by n the running index of the 
points separated by h. In correspondence with the lattice variable n, we can introduce the real variables 
X = hn. 

We can define on the same lattice a set of slow varying variables by introducing a small parameter 
e = and requiring that 

Uj = e-'n. (6) 

This is équivalent to sampling points from the original variables which are situated at a distance of h 
between them and then setting them on a lattice of spacing h. The corresponding slowly varying real 
variables xj are related to the variable x by the équation xj = e'x. 

2.2 Expansion of slowly varying functions. 

Let us study the relation between functions living on the différent lattices defined in section 12.11 We 
consider a function f = fn defined on the points of a lattice of index n. Let us assume that /„ = 
5ni,n2,...,nic, i-s. / dcpcuds OU a finitc number K of slow varying lattice variables rij j — 1,2, . . . , K 
defined as in ©. We want to get explicit expressions for, say, fn+i in terms of gni,n2,...,nK evaluated on 
the points of the ni, n2, ■ ■ ■ , tik lattices. At first let us consider the case, studied in when we have 
only two différent lattices, i.e. K = 1. Using the results obtained in this case we will then consider the 
case corresponding to K = 2. The gênerai case will than be obvions. 

ï) K = 1 {fn = .gril)- In this case we use the foUowing resuit presented in 

A'= 5ni = ( , )9n.+^ = E -[ ^(*' ^) 

i=0 ^ ^ i=k 

Here the coefficients P{i, k) are given by 

i 

p{i,3)^Y.^^^tei, (8) 

where lû is the ratio of the incrément in the lattice of variable n with respect to that of variable ni . In 
this case, taking into account équation u) = N . The coefhcients Sf and (3^ are the Stirling numbers 
of the first and second kind respectively j2| ■ Formula Q allow us to express a différence of order k in the 
lattice of variable ni in terms of an infinité number of différences on the lattice of variable n. The resuit 
{Tj) can be inverted and we get: 

°° k^ 

^'•^"^E-[Q(*'^)^'5»i, (9) 

i—k 

where the coefHcients Q{i, k) are given by ||HJ) with lu = N^^ — e. 

To get from équations 10 and (jÏÏÏ a finite approximation of the variation of Çm = fn we need to 
truncate the expansion in the r.h.s. by requiring a slow varying condition for the function /„. Let us 
introduce the foUowing définition: 

Définition. The function fn is a slow varying function of order p iff A^^^ fn = 0. 

From Définition (2.1) it foUows that a slow varying function of order p is a polynomial of degree p in 
n. From équations Q and © we see that also the foUowing statement holds: 

Theorem. /„ is a slow varying function of order p iff A^"*"^ g^ = 0, namely g^ is of order p. 

Equation 10 provide us with the formulae for fn+i in terms of g^ and its neighboring points in the 
case of slow varying functions of any order. Let us write down explicitly thèse expressions in the case of 
gni of order 1,2 and 3. 

9 p = 1. The formula © reduces to 

Afn = TfA^m, 
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i.e. fn+i reads 
p = 2. From équation we get 



fn+l = 5ni + -^{gni + l " ffm) + 0(iV 2). 



1a 1 a2 



and thus fn+i reads 

fn+i = 5ni + ^(-g«i+2 + 4g„i+i - 3.g„J + ^^(g«i+2 " 2.g„i+i + g„J + 0{N ). (10) 

• p = 3. From équation ^ we get 

1a 1-^a2 {1-N){1-2N) ^3 

and thus fn+i reads 

fn+i = ffm + g^(2gni+3 - 95„i+2 + 13.g„i+i - 6.g„J + 

1 / N 

+ ^^(-.9ni+3 + 45„i+2 - 55„i+i + 2g„J + 

+ ^^(3ni+3 - 3.grii+2 + SÇn. + l - .gnj + 0{N^^). 

In the next sections we will consider mainly the réduction of integrable discrète équations and we 
will be interested in obtaining from them integrable discrète équations. It is known |13j that a scalar 
differential-difïerence équation can possess higher conservation laws and thus be integrable only if it 
dépends symmetrically on the discrète variable, i.e. if the discrète équation is invariant with respect to 
the inversion of the lattice index. The results contained in do not provide us with symmetric formulas. 
To get symmetric formulas we start from équation Q and take into account the foUowing remarks: 

1. Formula ((TJ holds also if h is négative; 

2. For a slow varying function of order p, we have A^fn — A^fn+i, for ail £ E 7^. 

When /„ is a slow varying function of odd order we are not able to construct completely symmetric 
derivatives using just an odd nuniber of points centered around the ni point and thus /„±i can never be 
expressed in a symmetric form. 

Using the above remarks we can construct the symmetric version of (|10|) . From {Tj) we get: 

5„,+i = 9n, + NAf„ + In{N - l)A2/«, (11) 
where A^/n = fn+i — 2f„ + fn-i thanks to the remark [3 Using the remark QJwe can also write: 

9m-i - 9n, + iVA_i/„ + ^N{N - l)AVn, (12) 
where A_i /„ = /„_i — /„. From équations l(TÏ|) and ifT^ we obtain the foUowing form for fn+i- 

11 

fn+l = gm + ^(ffni+l - Sni-l) + ^^(.9ni+l " 2gn, + 9n^-l) + 0{N ). (13) 



II) K = 2 {fn = 9ni,n2)- The dérivation of the formulae in this case is done in the same spirit as for 
the symmetric expansion presented above, see équation (|13l) . Let us just consider the case when p — 2, 
as this is the lowest value of p for which wc can consider /„ as a function of the two scales ni and «2- 
From équation Hll|l we get: 

gn, + l,n, = 5ni,»2 + ^Ai/„,„ + ^N{N - l)Alfn,n, (14) 
gnun2 + l = 9n,.n2 + ^2În^n + ]^N^ [N^ ~ l)A^/„,„. (15) 
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Here the symbols Ai and A2 dénote différence operators which acts on the first and respectively on the 
second index of the function /„^„ = 5ni,n2: e.g. Ai/„,„ = /„+i,„ - /„,„ and A2/„,n = fn.n+i - în,n- 

Let us now consider a function gni,n2 where one shifts both indices by 1. From équation H14|l . taking 
into account that, from équation lO, for example, gni+i,n2 = fn+N,n, one bas: 



ffni + l,n2 + l — 9ni,n2 + l + ^^l/ri^n+W^ + -^N{N - l)Af/„ „+jv2, 

and using the resuit IjlSf) we can write équation (|16|l as 

1 



.9ni + l,n2 + l 



(16) 



(17) 



-NAi 



/„,„ + N^A2fn,n + ^N^iN^ - l)Alfn,n 



+ -NiN~l)Ai 



/„,„ + N^A^fn^n + -N%N^ - l)AlU 



g„,,„, + N'A2fn,n + -N\N' - l)A^/„,„ + 



+iVAi/„,„ + N^A^A^fn^n + -N^N^ - 1)AiA2/„,„ + 

+ ^N{N - l)Alfn,n + iV3(7V - l)A?A2/„,„ + ^N^'iN^ - 1){N - 1)A2a2/„,„. 

As, using the second remark, the second différence of dépends just on its nearest neighboring points, 
the right hand side of équation dépends, apart from = gni,n2^ on /„,„+!, fn,n-i, /n+i,n, /n-i,n, 
/n+i,n+i, /n+i,n-i, /n-i,n+i, and /„-i,„-i, 1-6. 8 unknowns. Starting from équations (O, and 
()17|l we can write down 8 équations, using the first remark, which define gm+i.ng' 5ni-i.n2 7 5ni,n2+ii 
5ni,n2-ii .9ni+i,«2+i> ffm+i.ns-i) S'ni-i,n2+i: and g„i_i,„2_i in terms of the functions fn+i,n+] with 
(i, j) — 0, ±1. Inverting this System of équations we get /„±i in term of gm,n2 and its shifted values: 

/n±l ~ 3Tii,ri2 i ^^(â'ni + l,n2 ~ .9ni-l,n2) + 27V^ (5ni + l,n2 ~ '^gni,n2 +5ïii — l,n2) + 



± 



+ 



1 



27V2 
1 

4ÎV3 



(5ni,«2 + l 5'ni,"2 — l) ^" 

(5ni+l,Tî2 + l ~ 5'ni~l,n2 + l " 5r!i + l,n2-l + .9ni-l,n2-l) + 0{N "^). 



(18) 



It is worthwhile to notice that the two lowest order (in N~^) terms of the expansion (|18() are just the 
sum of the first symmetric différences of gm and gn^. Thus in the continuons limit, when we divide by h 
and send h to zéro in such a way that x = hn, xi = hni and X2 — hn2 be finite, we will have 

f,x = egxi + e^gx2- 

Extra terms appear at the order and contain shifts in both ni and n2. 

When /„ is a slow varying function of order 2 in ni it can also be of order 1 m n2- In such a case 
équation p5|) is given by 



2 + 1 — 



A^A2/„ 



(19) 



Starting from équations p4|) , p9|l and a modified (|17|l we can get a set of 8 équations which allows us to 
get fn±i in terms of gm,n2 and its shifted values. In such a case /„±i reads 



/n±l ~ 5ni,n2 ^ ^^"i +-'-'"2 5iîi-l,n2)^' j^2^3ni ,n2±l 5ni,n2)"l" 



(20) 



It is possible to introduce two parameters in the définition of ni, 712 in terms of n. Let us define 



ni 



nMi 



n2 



nM2 
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where Mi and M2 are divisors of N and TV^ so that ni and n2 are integers numbers. In such a case 
équation lfTH|) reads 

_ I ^1 / \ ( 

Jn±l — 9ni,n2 ^ '2^\9ni + l,n2 ~ ffni - l,n2 ) + 2A^2 (5ni + l,n2 ~ ^ffni ,«2 "I" 5"! — l,n2 / "1" 

H ^^y3^(3ni + l,n2 + l ~ 5ni-l,n2 + l - 5ni + l,n2-l + ffni-l,n2-l) + ^(iV (21) 

and équation (|20|) accordingly. 

When we consider partial différence équations we have more than one independent variable. Let us 
consider the case of two independent lattices and a function fn,m defined on them. As the two lattices 
are independent the formulae presented above apply independently on each of the lattice variables. So, 
for instance, the variation /n+i,m when the function fn,m is a slowly varying function of order 2 of a 
lattice variable ni reads 

fn+l,m — Onum + ^^(Sm+l.m — 5ni-l,m) + ^^(5ni + l,m " "^grium + ffni-l,m) + 0{N ^). (22) 

A slightly less obvions situation appears when we consider /n+i,m+i, as new ternis will appear. We 
consider here just the case we will need later when 

Min M2TO n 

«1 = ^, = ^2 = j^. (23) 

If /ri. m is a slow varying function of first order in TO2 and of second order in both ni and mi, from 
équations (|21|l and (|22|) the variation fn+i,m+i reads 

/n+l,m+l — gni,mi,m2 + (ffni + l,mi ,m2 ~ .Çrii — l.mi ,r?j.2 ) + (.Çni ,mi + l.m2 ~ 5"! ."H - 1:™2 ) ~^ 

"I" ^^"i+l-^i'^^ .9rii — l,mi,m2 ^9ni,mi,m2) 

"I" 2^2 ^^"1 + 5"! ,mi — l,"i2 2f;„j^mi ,?Ti2 ) "1" 

H ^JY2 (5ni + l,mi + l,m2 ~l" l,»"!— 1.™2 ~ 91^,1+1,^^,1—1,^^,2 ~ .9ni - l,mi + l,m2 j + 

+ ^(5ni,mi,m2 + l ~ 3ni,mi,m2) + 0(A^^^). 



3 Multiscale réduction of nonlinear partial différence équations 

In the foUowing we will apply the formulae obtained in section 2 to some well known partial différence 
équations. Some of those are known to have a Lax pair and are associated to integrable partial differential 
équations. Others are concocted so as to have a real dispersion relation but with no particular reason 
why they should be integrable. The integrable équations we will consider here, the lattice modified 
KdV (mKdV), presented in section IXÏl the Hietarinta équation, presented in section and the lattice 
Volterra-Kac-Van Moerbeke (VKVM) équation, presented in section 13.31 are defined on four lattice 
points and are PAEs consistent around a cube 5 . From this property one can dérive their Lax équation. 

The lattice mKdV is an integrable équation of the same class of the lattice potential KdV and KdV 
[3] and it possesses a Lax pair ^U] . As from KdV we get by multiscale réduction the NLS j4: , the same 
we may expect here. To get an integrable discrète équation we expect a resulting discrète équation which 
is somehow symmetric. At least when ht with t ~ mht the differential différence équation we obtain 
must be symmetric in terms of the inversion of nj, i.e. if it contains rij+fe it will contain also nj^k- 

The non integrable KdV équation presented in section |^31is obtained by a straightforward discretiza- 
tion using the symmetric représentation of the derivatives so as to get a real dispersion relation. 

In ail the cases considered we will expand the solution of the nonlinear lattice équation around a wave 
solution of the linear part. In doing so we require that the wave solution be always bounded so that a 
perturbative expansion with slowly variable coefficients is meaningful. This can be always achieved if the 
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dispersion relation is real. This is always true if the équation can be rewritten in ternis of symmetric 
derivatives. Moreover to get a meaningful réduction we need to have a non trivial nonlinear dispersion 
relation. 



3.1 Réduction of the lattice mKdV 



The discrète analogue of the modified Korteweg-de Vries (mKdV) équation is given by the foUowing 
nonlinear PAE ^ i 



(24) 



This équation involves just four points which lay on two orthogonal infinité lattices and are the vertices of 
an elementary square. In équation H24I) Un^m is the dynamical field (real) variable at site (m, n) G Z x Z 
and p, g G R are the lattice parameters. Thèse are assumed différent from zéro and will go to zéro in the 
continuons limit so as to get the continuons mKdV. 

Carrying out the change of variable it„_m ^ 1 + Un,m, one can separate the linear and nonlinear parts 
of équation (|24|l : 

P {Un,m Wn,m+1 ^n+l,ï7i ^n+l,m+l) Q (^n.ïri 4" ^ïi+l,m ^n,m+l '^n+l,m+l) — 
= Q {Un^m Un+l.m ^ ''^ri,m+l ^n+l.m+l) ~ P (^n,m î^n,m+l ~ '^n+l,m ^n+l,m+l)- (25) 

Let us consider the linear part of équation H25|) . namely 

P {Un,m + Un,m+1 — Un+l,m — Un+l,m+l) — Q {Un,m + Wn+l,m — Un,m+1 — Un+l,m+l) = 0. (26) 

Given any initial condition u„.o the gênerai solution of équation H26|l is given by 

'{p-q) - {p + q)z' 



~ 27ri ^ 



(27) 



(p - q)z -{p + q) 

j— ' — ' 

Equation H27(l can be rewritten in a more natural way (from the continuons point of view) by defining 

{p-q)-{p + q)z 



{p - q)z -{p + q)' 

In such a case the solution l|27|l is written as a superposition of linear waves 

p i[kn—u{k)m\ n Qin 

The dispersion relation for thèse linear waves is given by 



Lj = —2 arctan 



- tan i — 
2 



(28) 



(29) 



(30) 



the same as for the lattice potential KdV (pKdV) équation 
differentiation with respect to k, we get the group velocity uj^k- 



OJ,k 



Apqz 



From équations (|28|l and H30() . by 
2pq 



[{p - q)z - {p + q)] [{p + q)z - {p — q)] p'^ + q"^ — (jP' — q^) cos k 



(31) 



The linear part of the PAE H25|l . i.e. équation H26|l . is solved in ternis of harmonies (|29|l if oj is given 
by H3()|l . The nonlinearity will couple the harmonies. This suggests to look for solutions of the PAE (j25|l 
written as a combination of modulated waves: 

oo oo 

Un,ra = ^ e^^i^^^ [En^mY + Y. '"'^nl (^n,™)^ (32) 
s=0 s=l 

where the functions ipn^m are slowly varying functions on the lattice, i.e. ipn)m = ipni,mi,m2 and — N~^. 
By b we mean the complex conjugate of a complex quantity b so that, for example, -Ë„,m = (-En,™)"^- 
The positive numbers /3s are to be determined in such a way that : 
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1- /?! < /3s V s = 0, 2, 3, . . . , oo. In gênerai it is possible to set (3i = 1. 

2. In the équation for ijjn^m = ijM,m, we require that the lowest order nonlinear ternis should match the 
slow time derivative of the linear part after having solved ail linear équations. This will provide a 
relation between 7 and the f3s. 

The fact that the second summation in équation H32I) starts from 5 = 1 and contains the complex 
conjugates of the terms of the first summation is due to the reality condition for the solutions of the P AE 

After introducing the expansion H32() in the PAE (|25|1 and analizing the coefficients of the varions 
harmonies {En^mY for s = 1, s = 2 and s = (as, assuming that /3s increases with s, the nonlinear terms 
will dépend only on the lowest s terms) we came to the conclusion that we can choose 

7 = 1, /3o = 2, ps = s, s>l. (33) 

The discrète slow varying variables ni, mi and m2 are defincd in terms of n and m by équation (|23|l . 

Having fixed the constants (3s according to équation (|33|l we can introduce the ansatz (|32|l into équation 
()25|l and get the determining équations. 

For s = 1 we get, at lowest order in e, 

V'ni.mi^msK'? --P)(l - Z n) ~ {p + q)(n - z)] ^ 0, 

which is identically solved by the dispersion relation H28|l . 
At we get the linear équation 

'0ni — l,mi ,m2 ) ~^ 

+ M2^[(j} - q)z - + q)]{^pn^^„i^ + i^„i2 - ^7ii,nn-l,m2) = (34) 

whose solution is given by 

provided that the integers Mi and M2 are choosen as 

Mi^Sn[ip~q)z-ip + q)i M2^Sz[ip~q)n + ip + q)], (35) 

where S" G C is a constant; S cannot be completely arbitrary since Mi and M2 are to be integer numbers. 
We will show in Appcndix A how it is possible to choose the complex constant S in such a way that Mi 
and M2 are in fact integer numbers as required by équations (|23|l . Substituting the expression of Çl given 
in H28|) into équation H35|l . we can rewrite Mi and M2 as: 

Mi^~S[{p + q)z-{p-q)], M2^j-—^^^ -y (36) 

[{p + q) - z{p - q)\ 



From équations (|31|l and H36|l we get: 



= ^, (37) 



i.e. the ratio M2/M1 is the group velocity. As Mi and M2 are integers, it foUows that not ail values of k 
are admissible as uj^k G Q- Let us notice that also 712 = «i + mi solves équation l|34(l by an appropriate 
choice of Mi and M2 ■ 



At we get a nonlinear équation for 4>n2,m2 which dépends on ^l^^ 



2, "12 • 



^n2,m2 + l Ç'n2,m2 ^1 (Ç^n2+2,m2 ^~ Ç^n2— 2,m2 ^ ^n2,m2 



+ C2 

where 



«2 + 1. "î2 ^ 4'n2 — l,m2 2 0„2,m2) ^5 '4'n2 ,in2 4'n2,m2 — 0) (38) 

/ C.2 2 (p - ç) - (p + q)z 

[{p - q)z - {p + q)] 

ry , .^2 {p + q){l + z^)-2{p-q)z 
C2 = 2pq{p-q)S z ^ -2 , (39) 

[{p-q)z- [p + q)\ 



C3 



2pq (p2 — ç2~j — z'^Y' 



z [{p ~ q)z - {p + q)]^ [ip + q)z -{p- q)f 
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Using the form of the complex constant S obtained in Appendix A, the coefRcients ()39fl read 

Cl = - ^\ 1^ — — [{p + q){cosk + ismk)-{p-q)], 
Ibpq 

C2 = ^^^^^[{p + q)c08k-{p-q)], (40) 
2pq{p^ — q^) sva^ k 

C3 ^ ' 



+ q'^) ~ {p^ — q"^) cosfc]2 

The coefficients H4()|l dépend on the integer constant M2. The integer Mi is then written ont in terms of 
M2 and reads 

2pq 

so that not aU values of k are admissible as Mi must be also an integer. See Appendix A for détails. 
The lowest order équations for the harmonie s = 2 appear at and give 

é^"^^ --(é ? 

It is easy to see that the choice implies that the coefficients of ail other harmonies are expressed in 
terms of 0„2,m2 and 0n2,m2- 

Taking thèse results into account the nonlinear équation PAE H38|l for 4>n2,m2 reads: 



i('/'n2,m2 + l ~ 0Ti2,m2) — C*! (0n2+2,m2 + 4'n2-2,m2 ~ 2 0ri2,m2) + 

«2+1, m2 4'n2 — l.m2 2 0„2,m2 «2, "12 I0n2,™2 r, (41) 

where Ci = —ici, z = 1, 2, C3 = —103/2, and the coefficients q's are given by équation (|40|) . It is easy to 
see that C3 is a real coefficient. 

The PAE (|41() is a completely discrète and local NLS équation depending on the first and second 
neighboring lattice points. At différence from the Ablowitz and Ladik ^ discrète NLS, the nonlinear 
term in H41|) is completely local. The PAE H41|l has a natural semi-continuous limit when m2 — > 00 as 
H2 — > in such a way that t2 = m2 H2 G M is finite. Setting n2 = n and t2 ^ t one gets the following 
nonlinear differential-difference équation: 

i = Cl {4>n+2 + 4>n-2 - 2 0„) + 6*2 (0,1+1 + 0,1-1 - 2 0,,) + C3 0„ |0np. (42) 

The continuons limit of the PAE H41() is obtained if we consider in équation H42|l the limit n ^ oo as 
_ffi ^ in such a way that x — nHi G K is finite. The resulting NLS équation reads 

i0,t = (4Ci+C2)0,.. + C3 0|0p, (43) 

where 

4 Cl + C2 = ^ ^ , ' . (44) 

4pq 

As the coefficient (|44f) is real, équation (|43ll is just the well known intcgrable NLS équation. 

3.2 Réduction of the Hietarinta équation 

In [S] Hietarinta introduces a new consistent around a cube PAE 

Un,m + 62 Un+l,m+l + «2 U,i+l,m + ^2 Mn,m+1 + «2 
Un,m + ei U,i+1,,„+1 + Oi Un+1,,„ + Oi Un.m + 1 + Cl 

where the four constants e^, G M, 1 < i < 2, are lattice parameters. 



= 0, (45) 
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By a direct calculation one can separate the linear and the nonlinear parts of the équation H45|l : 

Ol 02(61 ^2)Un.m + ei 62(01 - 02)Un+l.m+l + 
+ Cl 02(62 - Oi)u„+i,„i + 62 01(02 - 6i)u„,„i+i = 

~ [(02 ~ ei)u„+i^„i + (62 — Oi)u„^,„+i]u„^„i u„+i,,„+i + 

+ [{01 

+ [01 (62 - 02)Un,m + l + 02(01 - 6i)u„+i,m] U„,,„ + 

+ [62(61 - Oi)u„,„+i + 61(02 — e2)'^n-\-l,m] ■^n+l.m+l + 

H" (^2 62 Oi Cl)('îi7i^77T, 'Un-\-l,7n-\-l ^ '^n+l.m '^n,m-\-l 

). (46) 

Let us now solve the linear part of the PAE (|46|l : 

Ol 02(61 - 62)u,um + 61 62(01 - 02)u„+l,m+l + 
+ 61 01(62 - 02)lt„+l,m + 62 02(01 - 6l)u„,m+l = 0. (47) 



Defining 

^^gifc o J. ^-iu, _ 02 [ei(62 - oi)z + 01(61 - 62)] 

62 [61(02 - oi)z + 01(61 - 02)] ■ 
the (complex) dispersion relation for thèse linear waves is given by 

161 01(02 - 62) tan (fc/2) 



2 arctan 



[oi 02 (61 - 62) + 61 62 (oi - 02)] tan (fc/2) + ie2 02(61 - oi) 

The dispersion relation is a real function of k if the foUowing condition holds: 

Ol 02 (61 - 62) + 61 62 (oi - 02) = 0. (48) 

To give a meaning to the expansion (I32|l we have to require that the dispersion relation Lo{k) be a real 
function for fc real. From équation 148|) we can write 02 in terms of oi, 61, 62, and get: 

n = 6-- = (49) 

01(61 - 62)Z + 61(62 - Ol) 



so that the real dispersion relation reads 



2 arctan 



62(61 + Ol) - 2 61 Ol fk 
tan — 



62(01 - 61) 

From équation H49() . by differentiation with respect to fc, we get the real group velocity uj^k 



(50) 



^ ^ _ 62(01 - 6i)[2 6i Ol - 62(61 + Ol)] Z 

[61(01 - 62)2 + 01(62 - 6i)][0i(6i - 62)2; + 61(62 - Ol)] ■ 

The PAE H47|) lias a bounded wave solution given by équation (|29|l . where Vl is given by équation 
(|49|l . So we can look for solutions of the PAE (|46|l in the form of a combination of modulated waves 

. is) . . . . (s) (s) 

(j32|l . where the functions V'n,™ are slowly varying functions on the lattice, i.e. ipn'm = 4'ni,mi,m2 and 

Introducing the expansion H32|l in the Hietarinta équation (|46f) and considering the équations for 
s = 1, s = 2 and s = harmonies we deduce that the choice (|33|l is still valid. Moreover, the discrète 
slow varying variables ni , mi and 7712 are defined in terms of n and m by équation (|23|) . 

Having fixed the constants Ps we can now introduce the ansatz H32|l into équation H4t)|) and pick out 
the coefficients of the varions harmonies {En^mY to get the determining équations. 

For s — 1, having defined ipn}n = 'ipn,m, we obtain an équation at the first order in e which is identically 
solved by the dispersion relation H5()|l. 

At we get the linear équation 

Ml z [01(61 - 62)ri + 61(01 - 62)] (V'ni + l,mi,m2 " "^ni - 1,™! ,^2 ) + 
+ M2 [01(61 - 62)2 - 61(01 - 62)] (V'ni,mi + l,m2 " V""! ,mi -1,^2 ) = Oj 
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whose solution is given by 

provided that the integers Mi and M2 are choosen as 

Ml = SVl [oi(ei - 62)2 - ei(oi - 62)], A/2 = 5*2; [oi(ei - 62)^! + ei(oi - 62)], 



(52) 



where S* G C is a constant. Inscrting Vl given by équation (|49|1 in équation (|52l) we can show that the 
ratio M2/M1 coincides with the group velocity (|ïïT|l . 



At e'^ we get a nonhnear équation for (pn2,m2 which dépends on ïl^n2,m2 and ipn2,m2- It reads 



,,(2) 



'^«2 ,m2 + l 



12, "12 

n2,m2 ^1 (0n2+2,m2 4^n2~2,rn2 2Ç^n2,m2) 

2 I 



+ C2 {4^n2 + l,m2 Ç^n2 — l,m2 2(/>^2,m2) ^~ 0n2 ,m2 |Ç^n2 ,m2 
+ C4 '4'n},m2 ^n2,m2 + C5 1pn],m2 4'n2,m2 — 0; 



(53) 



where the coefficients 1 < i < 5, dépend on z, S* and the lattice parameters ei, 62, oi and are given in 
Appendix B as their expressions are rather compHcated. 

The functions il>n},m2 and '4>n},m2 that appear in équation (|53|l are obtained by considering the équa- 
tions for the harmonies s = 0, at the third order in e, and s = 2 at the second one. From them we get: 



-■Pl{(t>n2,m2f', 



rn2 + l,m2 



■^«2 — l,JTi2 



with 



^ P2 VPn2,m2\Çn2+l.m2 ' 
Cl Z — Oi 



Pl 



^«2— 1,™2) ~^ 4'n2,m2{4'n2 + l,m2 

ei + oi 



(54) 

^n2-l,m2)]: (55) 



ei oi{z - 1) 



P2 



ei oi 

/,(2) 



From équations (|54|) and H55|l we évince that both '0n2,m2 and '4'n2-m2 are expressed in term of 4>n2,m2- 
particular, we notice that ipn2,m2 dépends from 4'n2,m2 in a local way while ipn},m2 dépends from çi^^ 
in a non local way through a summation, namely 



In 



= (-1)" 
Yn2,m2 V / 



'«1 +P2^ (-l)''(0j\m2 0j + l,m2 + 0J,m2 '/'i+l.m2) 



J="2 



W2, 



(56) 



where wi,W2 are two arbitrary summation constants. 

Inserting 4'i^2,rn2 given by équation H54() in équation (|53|l we get 

Ç^n2,m2 + l 0n2,m2 ~^ '-^l (0n2+2,m2 ~t~ ^«2— 2,m2 ^ ^7i2,ra2^ ~t~ 

2(/'„2,m2) + C3ç!'„2:m2l'/'n2,m2p + C4 V'l°2m2 '?^"2:m2 = 0> (57) 



+ C2 ((pii2+l,m2 + <Pji2-l:m2 



where, using the form of the complex constant S, see Appendix A, and the fact that z = e''', the 
coefficients are: 



Cl = 
C2 

C3 

C4 



M. 



2 ; 



P2 [-Pl (cos A: + i sin A:) + P2] 

P2(PlCOsfc + P2) ,,2 

^^^^ ^' 

2(Pi - P2)[Pi(ei - 62) + P2(e2 - oi)](cosfc - 1) 

62 (Ol 62 + P2 ) (Pf + P| + 2 Pl P2 COS /c) 
2(Pl - P2)(Pl -6^+01 62)(C0S fc - 1) 



62(Pi2 +P| +2P1P2 cos/e) 



with 



Pl = 61(62 - Ol), P2 = 01(61 - 62). 

Here M2 is an arbitrary integer number, while Mi is given by (see Appendix A) 

P? +Pi + "2 Pl P2 COS k 



Ml = M2- 



p2 _ p2 
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3.3 Réduction of the lattice VKVM équation 

The completely discrète version of the Voherra~Kac-Van Moerbeke (VKVM) équation is given by the 
following PAE : 

= 7i \^°) 

Here a is a real lattice parameter and w„^m is a real field. 

The dispersion relation of the linear part of équation (|58|l is trivial. So, we carry out the change of 
variable Un,m 1 + u„.„i. Then one can split équation (|58|l into the linear and nonlinear parts: 

a(Un+l,m+l - Un,m) + (1 - a)(Wn+l,m ^ Un.m+l) = a(Un,m Uji+l,m - Wn+l,m+l Urî,m+l)- (59) 

The dispersion relation for the linear waves is given by 

(2a — 1) sin k 



a{z + l)-z 

" = — ; 7^ , = arctan ,„ „ „ , „ , n 

q;(z + 1) - 1 [ (2a^ - 2a + 1) cos A: + 2a(a - 1) 

From équations (|60|l . by differentiation with respect to k, we get the group velocity uoj^ 



(60) 



^ {2a-l)z (2a -1) 

[a(z + l)-z][a(z + !)-!] 2a(a - l)(cos fc + 1) + 1 ' ^ ' 

We now consider a solution of the PAE H59I) in the form of a combination of modulated waves, see 
équation l|32|) . where i?n,m is given by équation l|29|) with Çl as in équation H6UI) . As in the previous cases 

(s) . . . . (s) (s) _1 

the functions \l>n,m are to be slowly varying functions on the lattice, i.e. xpn'm — V'ni,mi,m2 s-nd = N . 

Introducing the expansion H32|) in the lattice VKVM équation (|59|l and considering the équations for 
s = 1, s = 2 and s = we deduce that the choice (I33|l is still valid. The discrète slow varying variables 
ni, nii and m2 are defined in ternis of n and m by the positions H23|) . where Mi, M2 £ Z. 

For s = 1 we obtain an équation at the first order in e which is identically solved by the dispersion 
relation H60() . 

At we get a linear équation 

Ml Z {a[VL - 1) + 1] (V^„i + l,mi,m2 - -iAni-l^mi.ma) + 
+ M2 fi [a(z + 1) - 1] (■(/'ni,mi + l,m2 " V'rii ,mi -l.ma ) = 0, 

whose solution is given by 

^^711,7711,7712 ^772,77725 ^1 ?Til. 

provided that the integers Mi and Af2 are choosen as 

Ml = 5rî[a(z + 1) - 1], Af2 = S'z[a(r2-l) + l], (62) 

where S* G C is a constant. Inscrting $7 given by équation H6U|I in équation (|62|) we get that the ratio 
Mil Ml coïncides with the group velocity H61|l . As shown in Appendix A it is possible to choose the 
complex constant S in such a way that Mi and M2 are in fact integer numbers. 

At e" we get a nonlinear équation for (j)n2,m2 which dépends on V'n2'',77i2 and 'ip'n2,m2- It reads: 

'/>772, 7772 + 1 " ^772 , 77l2 + ^1 ( (/>772 +2 , 7772 + </'772 - 2 , 7772 ^ ^ (/l,i2,77l2) + 
+ C2 (0772 + 1,7)72 + 0772-1,7712 ~ ^ 4>n2,m2) ~l" ''S '4''^n},m2 4'n2,m2 + C4 'ipl^2,m2 <^772,7r72 = 0, (63) 

where the coefficients c, , 1 < î < 4, are: 



Cl = a 5 z 



2 2 (l-2a)[a(z + l)-z] 



4[a(z + 1) - 1]2 



2 (2a-l)[a(z + l)2-z2_i] 
C2 = a D z — 



C3 



C4 



2[a(z + l)-l]2 
a(l-z2) 



[a(z + l)-z][a(z + !)-!]' 

a(l-z2)(z2-z+l) 
[a(z + 1) - z][a(z + 1) - l]z' 
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The functions ipn^^m^ and ^''n},m2 that appear in équation are obtained by considering the équa- 
tions for the harmonies s = 0, at the second order in e, and s = 2, at the third one. We get the following 
équations: 

V-i^Us =Pl('^"2,m2)^ (64) 
V'ia + l.ma ~ '^n2-l,m2 ^ ^2 [^na ,m2 ('/'na +I,m2 ~ 0n2~l,m2) + '/'n2,ni2 (^"a + l^ma " ^«2-1,^2)]: (65) 

with 

_ (l-2a)z _ 2a(l + z2) - z2 - 1 

" a(z + l)2-z2-l' " (l-z2)(a-l) ■ 

From équations (|64() and H65|l we évince that both ^n},m2 and '0i2'',"i2 are expressed in term of 4>n2.m2- In 
particular '4'n}.,m2 admits a non local expansion as in équation (j56|l . Inserting ^/'i2^,m2 given by équation 
(|64|l in équation Ht)3|) we get 

0n2,m2 + l 0n2,m2 "t" (^^712 H-2,m2 "t" 0ïi2— 2,m2 4^n2.m2^ ~^ 

+ C2 ((/'„2 + l,m2 + 0«2-l,m2 - 2 0„2,m2) + C3 ^'^"^ma '/'n2,m2 + C4 ^«2 ,m2 [(/-na ,m2 1^ = 0; (66) 

where C4 = c^pi- 

Let us Write the coefficients that appear in équation 166|l . i.e. ci, 02,03,64, using the form of the 
complex constant S, see Appendix A, and the fact that z ^ . We get 



(a — l)(cos k + isin k) -\- a 
4(2a- 1) 



Cl = —a jYP. 7^ -^^2 I 



(a — 1) cos + a , 
2a-l ^ 
2a sin fc 

C3 = -i 



C4 = 1 ■ 



2a(a - l)(cosfc + 1) + l' 

(2cosfc — 1) sin/c 



[a - 1) (cos A; - 1) [2a(a - 1) (cos fc + 1) + 1] ' 
As in the previous cases we can choose the integer number M2, while Mi is given by (see Appendix A) 

2a{a - l)(cosÂ; + 1) + 1 



Ml = M2- 



2Q! - 1 



3.4 Réduction of a non integrable lattice KdV équation 

Let us now consider the following non integrable lattice KdV équation: 

Un,m+1 — Un,m-1 — ^(Un+3,m " 3u„+i^„i + 3u„_i^„i — U„-3^m) + /3[(u„+l,m)^ " {Un~l,mf']i (67) 

where a, /? G K are the lattice parameters and u„.m is a real field. 

As we did for previous cases we apply the standard discrète Fourier transform procédure introducing 
Un,m — 2" î^™ into the linear part of the lattice équation H67I) . Here z = and fi = e~"^. We easily 
get: 

Hence the dispersion relation reads 
and the corresponding group velocity is 



Q. - Q.-^ ^ j{z ~ z'^f . (68) 



uj — arcsin (a sin'^ fc) 



3 aH, (z* — 1)(22 — 1) 3acosfcsin^fc 
"■^■ = -4TW P = Vl-a^sin^- ^ ^ 

Introducing the expansion l|32|) into the PAE ||ïï7Jl, where En^m is given by équation H29|l and taking 
into account that 

k k^ 

fn±k = ± -^{gni + l - grn-l) + |^(5ni + l " 2.g„i ) + 0(iV"^), 
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we get the standard choice (|33f) . 

Let us consider now the équations for the harmonies s = 1. The équation at the order e is identieally 
satisfied by taking into account the dispersion relation (|68|l . The équation at the order is satisfied if 
we introduce the index n2 = ni — toi when Mi and M2 are choosen as 

«, = 5(n + i), ...,-^5. '-'-'y-" . (70) 

We notice that the group velocity ui^k (I69|l coincides again with the ratio M2/M1, as in équation 1)37(1 . 
From équation H7U|) . using the fact that z = e'*^ and fl = e^"^ , wc obtain 

Ml = — 25cosijj, M2 = — 65*0? cos fc sin^ fc. 

We can now fix the (real) constant S* in a such a way that Mi is an integer number; M2 wih be an integer 
if the group velocity is a rational number. Hence not ail values of k are admissible, but only those which 
make uj^k l|69|l rational. 

The équation at the order e'^ is given by 



+ Cl (0„2+l.m2 + '/'ri2-l.m2 ~ 2(/'„2.m2) + C2{(f>n2,m2 (f'nlm2 + '^n2,ni2 0n2^,m2 ) = 0, (71) 



where 4'n2,m2 = '0ni,mi,m2i '^2 = ni — Toi and ci,C2 are known, easy to compute but too complicate to 
Write down, complex coefficients depending on z and on the lattice parameter a. The functions 4'n2,m2 
and i}jn2,m2 that appear in équation H71|l are obtained by considering the équations for the harmonies 
s = 0, at the third order in e, and s = 2, at the second one. We get the foUowing équations: 

^i'U =Pi('^»2,™2)', (72) 

''/'n2 + l,m2 ~ ''J^n2-l,m2 ~ P2 [0^2 ,m2 ('?^n2 + l,m2 ~ ^«2 - l,m2 ) + <^n2 ,ni2 (</'n2 + l,m2 ~ 4'n2-l,m2 )], (73) 

where pi , p2 are known complex coefficients depending on z and on the lattice parameters. From équations 
(|7^ and l(7ÏÏ| we évince that both 4'n},m2 and Tpn2,m2 are expressed in term of (/)„2,m2- in the previous 
cases 1^122,^12 admits a non local expansion as in équation (|56|l . Hence the PAE (|71|l is a well defined 
lattice équation in the field variable 4>n2,m2- 



4 Conclusive remarks 

In this paper we have shown that we can construct a well defined procédure to carry out the reductive 
perturbation technique on the lattice. In this case, at différence with respect to the differential-difference 
case, we are able to solvc ail linear équations and thus can obtain a final nonlinear différence équation. 
To do so we had to apply some non trivial but at the end obvions tricks which consist in the introduction 
of appropriate lattice variables so as to be able to perform the symmetric réduction of the linear discrète 
wave équation. 

Applying the perturbative-reductive technique to some integrable and non integrable équations we 
obtain some new completely discrète NLS équations. As some of thèse équations H41|l . H57|l and (|66|l 
come from the réduction of integrable équations we expect them to be also integrable. However they are 
very différent from the Ablowitz-Ladik discrète-discrète NLS ^ as ail contains, apart from the nearest 
neighboring points, also the points n ± 2 and either they are completely local or they have non local 
completely irregular terms (depending on (—1)"). 

So we are at the moment, from one side extending our analysis to other well known integrable équa- 
tions, like the discrète time Toda lattice, the sine-Gordon and the Volterra équations and from the other 
using the integrability properties of the starting nonlinear équations (i.e. Lax pairs or generalized sym- 
metries) to show the integrability of the derived équations. If our équation are integrable than we have 
presented a very important tool for obtaining new integrable équations and for analyzing the far field 
behavior of physical problems described by differential-difference or partial différence équations. 

In the dérivation we introduced the request that the far field expansion of a slow varying function on 
the lattice should dépend on the discrète asymptotic variables in a symmetric way. As a conséquence of 
this ansatz we got that the non local resulting équation dépends on (—1)". This may not be a necessary 
ansatz and work is in progress in this direction. 
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Appendix A 

From équations H35|) . H52I) and (|62|) we get that the coefficients Mi and M2 can always be written in the 
foUowing gênerai form: 

Mi^ Sn{Pz-Q), M2^Sz{Pn + Q), (74) 
where 5 G C is a suitable constant such that Afi, M2 G Z, il = e^"^, z = e'*^ and P,Q G R are given by: 

P ~ p — q, Q = P + q lattice mKdV équation, 

P — oi{ei — 62) , Q = ei(oi — 62) Hietarinta équation, 

P — a, Q — 1 — a lattice VKVM équation. 

The (real) dispersion relation for the linear parts of the above lattice équations can be written in the 
form 

Let us now define the complex constant S" as S* = p e'^, with p G and — tt < 6 < tt. From équations 
(dj) and ^ we get 

Re(Mi) ^ p[P cos{9) ~ Q cos{e + k)], (76) 
Im(Mi) =p[Psin(6') -Qsin(6' + fc)], (77) 

- P^+y-'2PQcosfc [^--(^) - + ^''^ 

Since Mi,M2 G Z we have to require that Im(A/i) = Im(A/2) = 0. From équations H77|) and (|79|l we 
obtain 

/ Q gjjj h \ 

6> ^ - arctan „ , ^ +^71", ^ G Z. (80) 



] cos k — p ^ 

We have now to require that Mi = Re(Mi),M2 = Re(Af2) G Z. According to équations (|75|l . (|78|l and 
(pHl we get 

Ml = {-\fp{P^ + Q^-2PQcoskf'^, (81) 

^^^^(-i)^^(P^T^^2î|W^- («2) 

We can fix arbitrarily the integer number Af2, (or equivalently Afi) and express Mi (or M2) in ternis of 

From (|83|l we can see that not ail values of k are admissible since Afi has to be integer. 

Let us finally notice that the fact that équation (|83|l contains cos(fc) implies that the ratios M2/M1 
and P/Q are constrained as —1 < cos(A;) < 1. We have the foUowing cases (see Figure 1): 

P/Qg(1,«3) ^ A/2/A/1 G [(P/Q- 1)/(P/Q + 1), (P/Q + l)/(P/g-l)] G Q, 

P/Q G (0,1) ^ Af2/Afi G [(P/Q + 1)/(P/Q-1), (P/Q- 1)/(P/Q + 1)] G Q, 

P/Q G (-1,0) => M2/Afi G [(P/Q-1)/(P/Q + 1),(P/Q+1)/(P/Q-1)] gQ, 

P/Q G (-00,-1) => A'/2/Afi G [(P/Q + 1)/(P/Q-1), (P/Q- 1)/(P/Q + 1)] G Q. 
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A Il/M, 




P/Û 



Figure 1: The grey zones dénote the allowed régions for the ratio M2/M1 in ternis of the ratio P/Q 



Appendix B 

The coefficients 1 < î < 5 that appear in équation (|53|l are 

ç2 P2(Pf-P|)[A(l + ^')+2P2^] 



C3 



C4 = 



C5 



2(Pi + P2zy 

(z ~ l)(Pi - P2)[Qi^^ + Q2Z^ + Qsz^ + Q4^' + Q5^ + Qe] 

e2{Pi-eie2)iPi+P2z)HPiz + P2)^z 
(z-l)^(Fi-P2)(ei-ei 62 + ^2) 

62(^22 + Pl)(PlZ + P2) 

(z - l)^(Pi - P2)[Riz* + ^2^^^ + fla^ + Ri] 
e2{Pi+P2zf{Piz + P2Yz 
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with 

P\ = ei(e2 - oi), 

< 

Pi = oi(ei - 62), 

v 

Qi =PiP2(Piei + P2e2), 

Q2 = Pf (ei - es) + P|(e2 - oi) + Pi P2(P2 ei + 2Pi - Pi oi), 
Q3 = -A[A'(ei - 62) + P|(ei + 4oi - 3 63) + Pi P2(3e2 - ei)], 

< 

Q4 = -P2[p2(4ei - 362 + 01) + P|(oi - 62) + Pi P2(3e2 - oi)], 
Qs = -Pi (ei - 62) - P|(e2 - 01) - Pi P2(P2 ei - 2Pi 62 - Pi d), 

' Ri = P2[P? + Pi + P1P2 + P2{el - ei 62)], 
P2 = Pi + (ei - ei e2)(P| - Pf ) + Pi P2(ei - ei 62 + Pi + 3P2), 

< 

iÎ3 = -Pi - (ei - ei e2)(P| - P^) - Pi P2(ei - d 62 + Pi + P2), 
^ P4 = Pi(Piei-P|-Pieie2). 
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